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1.  Modeling  Of  A  Piezoelectric  FGM  Plate 


Analytical  study  on  piezoelectric  laminate  with  functionally  graded  microstructure  (FGM) 
was  established  to  construct  a  hierarchical  modeling  (Almajid  et.  al.  (2000)).  The  study 
consist  of  two  models,  one  at  microscopic  level  focusing  on  the  microstructure  of  each 
FGM  layer,  the  other  on  laminating  FGM  layers  to  predict  the  out-of-plane  displacement 
of  the  FGM  plate  under  applied  electric  field.  The  standard  piezoelectric  actuator  plate  is 
of  bimorph  type.  These  kinds  of  actuators  develop  high  stresses  in  the  interface  between 
the  materials  as  shown  in  Figure  11.  An  FGM  model  has  been  investigated  to  reduce  the 
otherwise  high  stresses.  FGM  have  shown  a  remarkable  reduction  of  the  residual  stresses. 


Figure  1 1  Typical  stresses  generated  in  a  bimorph  piezo-actuator 


Figure  12  shows  typical  FGM  structure  where  the  mechanical  and  physical  properties  vary 
along  the  z-direction. 


z 


Figure  12  Piezoelectric  FGM  plate 
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A  piezoelectric  laminate  consists  of  n  laminae,  each  being  a  piezoelectric  material  with 
specified  electroelastic  properties.  The  constitutive  equations  of  a  piezoelectric  material  in 
the  absence  of  temperature  effects  are  given  by 

<7  =  C ■■„£,,  -e  ■  E  (4) 

wherecTy  are  the  stress  and  strain  tensor  components  respectively,  is  the  electric 

field  vector  component,  is  the  elastic  stiffness  tensor,  and  are  the  piezoelectric 

coefficients.  •  The  electroelastic  constants  of  each  lamina  may  be  computed  from  some 
micromechanical  model  such  as  the  Eshelby  type  method  as  outlined  by  Dunn  and  Taya 
(1993a,  b).  Two  analytical  modeling  has  been  established  to  predict  the  out-of-plane 
displacement  and  stress  field  of  laminated  piezoelectric  plate.  Classical  Lamination 
Theory  (CLT)  and  2D-elasticity  solution. 

1.1  Classical  Lamination  Theory  (CLT) 

Classical  lamination  theory  (CLT),  as  found  in  any  composite  textbook  such  as 
Gibson(1994)  assumes  a  state  of  plane  stress  along  the  z-axis  where,  a ^  =  0. 

The  constitutive  equations  of  a  piezoelectric  lamina,  eq.  4,  under  applied  electric  field  in 
the  z  direction  only,  and  under  the  assumption  of  plane  stress  along  the  z-axis,  reduces  to 
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It  is  noted  here  that  Q.j ,  e-j  are  the  reduced  stiffness  constants  and  reduced  piezoelectric 


constants  that  are  modified  by  the  assumption  of  plane  stress  and  where  £°,£°,  and 
are  the  in-plane  strain  components  at  mid-plane,  z=0  as  shown  in  Fig.  11,  x’^,x',,,and 
are  the  curvatures  of  the  plate. 
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The  resultant  in-plarie  forces  and  bending  moments  are  defined  by 


{W,M}=|;J{crX*.z*) 

'■='  '',-1 


(7) 


carrying  out  the  integration  through  the  plate  thickness  of  h,  the  resultant  forces  and 
bending  moments  can  then  be  written  as 


'n' 

’a 

B' 

K1 

B 

D 

lu’ 

where 

[A,B,D\^'^'^[Q\[dz,zdz,z'dz)  (9a) 

n  ^ 

[N,M  ]^  =  X  ^dz)  (9b) 

'='  A, -I 

where  /i,.y  and  hi  are  the  distance  from  the  bottom  of  the  plate  to  the  bottom  and  top 
interfaces  of  the  i'''  lamina  respectively,  with  i  =  1,2,  ...  n.  Hence  the  thickness  of  the  f'' 
lamina  is  then  hrhi.i.  Under  the  applied  electric  field  or  temperature  change  only,  i.e. 
M=N=0  in  eq.  (8),  the  in-plane  strain  and  curvature  can  be  solved  as 
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a  b 

K 
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(10) 


where 


a  b 
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b  d_ 

B 

D 

Under  a  given  electric  field  throughout  the  laminate,  one  can  then  predict  the  stress  and 
displacement  field  of  each  layer  as  well  as  the  out-of-plane  displacement  of  the  composite 
plate. 

Different  types  of  FGM  were  studied  as  shown  in  the  Figure  13  where  the 
conventional  bimorph  type,  one-sided  FGM,  and  FGM  bimorph  are  shown. 
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+v— 

(2)  bimorph  (b)  one  sided  FGM  laminate  (c)  FGM  bimorph 

Figure  13  Different  types  of  FGM 


The  distribution  of  mechanical  and  piezoelectric  (PZT)  properties  varies  throughout  the 
FGM  as  shown  in  Figure  14. 


z 


Figure  14  The  distribution  of  mechanical  and  piezoelectric  properties  in  the  FGM 


The  distribution  of  both  mechanical  and  piezoelectric  properties  could  be  explained 
through  the  following  formula 

C,=[0.l  +  0.225(i-l)]C,3, 

=[0.6  +  0.l(i-l)]e,2, 

where  C, ,  e-  are  the  mechanical  and  piezoelectric  properties  of  the  i-th  layer,  respectively. 
.C PZT  PZT  ^he  mechanical  and  piezoelectric  properties  of  the  piezoelectric  material. 
The  mechanical  and  piezoelectric  profile  of  the  FGM  bimorph  is  shown  in  Figure  15.  The 
properties  are  increasing  toward  the  middle  as  FGM-bimorph  type  A  while  they  are 
decreasing  in  FGM-bimorph  type  B. 
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The  mechanical  and  piezoelectric  properties  of  the  piezoelectric  and  the  mid-plane  gold 
layer  are  shown  in  Table  5. 


Cii 

(GPa) 

C,2 

(GPa) 

Cu 

(GPa) 

C3.2 

(GPa) 

C44 

(GPa) 

e.u 

(C/m-) 

^33 

(C/nr) 

ei5 . 

(C/nr) 

PZT 

121 

75.4 

75.2 

111 

21.1 

-5.4 

15.8 

12.3 

Au 

101.4.^ 

41.63 

41.63 

101.43 

29.9 

0 

0 

0 

Table  5  mechanical  and  piezoelectric  properties 


The  out  of  plane  displacement  of  each  case  were  calculated  and  tabulated  in  Table  6. 


Type  of  FGM 

microstructure 

Standard  bimorph, 

Fig.  2(a) 

FGM  one-sided, 

Fig.  2(b) 

Bimorph  FGM, 

Fig.  4(a)  Type  A 

Bimorph  FGM, 

Fig.  4(b)  Type  B 

Electric  Field 

(V/mm) 

1,000 

1,000 

1,000 

1,000 

Applied  Voltage  (V) 

208 

500 

208 

208 

Laminate  thickness 

(mm) 

0.5 

0.5 

0.5 

0.5 

No.  of  layers 

2 

5 

10 

10 

Layer  thickness(mm) 

. 

0.208 

Au  =  0.082 

0.1 

0.0416 

Au  =  0.082 

0.0416 

Au  =  0.082 

Curvature  (1/m) 

0.99 

0.992 

2.15 

0.773 

Max.  (MPa) 

1  12.18 

1  8.5 

9.8 

9.36 

Out-of  -plane 

displacement  (mm) 

for  plate  length  of 

50mm 

0.618 

0.62 

1.34 

0.48 

Table  6  Comparison  of  out-of-plane  displacement  and  maximum  in-plane  stress  for  various  types  of  piezo¬ 
composite  plates 


The  FGM  bimorph  with  mid-plane  layer  exhibited  the  maximum  out  of  plane 
displacement.  The  results  of  the  in-plane  stresses  between  layers  of  each  case  have  been 
are  shown  in  Figure  16. 
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z/h 

(a)  FGM  bimorph  type  B 


z/h 

(b)  FGM  bimorph  type  A 


.(MPa) 


z/h 

(c)  rainbow  type 


z/h 

(d)  standard  bimorph 


Figure  16  in-plane  stresses  throughout  the  plate  thickness  (z/h) 


As  it  can  be  seen  from  the  graph,  there  exist  high  in-plane  stresses  between  the  layers. 
That  stress  reached  high  in  the  case  of  bimorph.  The  high  residual  stresses  have  been 
reduced  through  the  use  of  FGM  especially  the  FGM  bimorph  with  intermediate  layer. 

The  gold  layer  serves  as  relaxation  for  the  piezoelectric  since  it  does  not  posses 
piezoelectric  properties.  There  is  a  need  to  investigate  on  the  shear  stresses  associated  with 
the  in-plane  stresses  by  using  a  more  accurate  model  such  as  FEM. 
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1.2  2D  elasticity  model  for  piezo  composite  laminates  in  cylindrical  bending 


z 


Figure  17  Simply  supported  piezoelectric  plate  under  a  state  of  plane  strain  in  the  y-direction 

For  a  laminate  of  multi  orthotropic  layers  of  piezoelectric  material,  the  cylindrical  bending 
was  constructed  by  following  the  solution  method  developed  by  Pagano  (1969)  and  using 
the  condition  of  plane  strain  along  the  y-direction  where,  =  0  as  shown  in 

Figure  17.  The  laminate  is  simply  supported  on  the  ends.  The  constitutive  equations  of  a 
piezoelectric  material  in  absence  of  temperature  are  given  by 

^ij  ~  ^ijkl^kl 

where (7, are  stress  and  strain  tensor  components,  E,„  is  electric  field  vector 
component,  Sy/.,  are  elastic  compliance,  df„„y  are  the  piezoelectric  coefficients.  The 
reduced  constitutive  equations  for  plane  strain  is  given  as. 


]  ^11  ]  ^  ^  ^31  p.v 

''6,  '  =  /?|3  /?33  0  <  <7,  ^ -h  0  0  z/33  <  Ey  >  (12) 

where  and  d^j  are  the  reduced  compliance  and  piezoelectric  coefficient  and  are 
defined  as 


The  boundary  conditions  of  the  simply  supported  piezoelectric  laminate  are  given  by 

a^{x  =  0,z)  =  (7,{x  =  l,z)  =  0 
w(jr  =  0,  z)  =  w{x  =  l,z)  =  0 
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while  the  prescribed  surface  traction-free  conditions  of  the  upper  and  lower  surfaces  of  the 
laminate  require 


The  interface  continuity  conditions  at  each  layer  are  given  by 


where  i  represent  the  interface  E)  number  between  the  f’  and  the  {i+lf'  laminae  with  i  = 
1  ~  N-1  and  the  topmost  and  bottommost  interfaces  are  defined  as  i”  and  (N-lf' 
interfaces,  respectively,  while  hi  represent  the  thickness  of  the  /'*  laminae.  The  top  layer 
corresponds  to  i  =  1  while  the  bottom  layer  is  i=N.  The  electric  field  is  applied  in  the  z- 
direction  through  the  thickness  of  the  laminate.  The  electric  field  is  expressed  in  the  form 
of  Fourier  series  as 


E.  =  E^,  —sin (Ax) 


E.=E,.=0 


(17) 


where  A  =  — ,  and  for  higher  values  of  n,  a  nearly  constant  electric  field  can  be  realized. 
Following  Pagano  (1969),  the  stress  field  in  the  i"'  layer  is  assumed  to  be 
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where  f\  (z)  is  an  unknown  function  to  be  determined.  The  equilibrium  equations  given  by 
=0  (19) 

are  automatically  satisfied  if  eq.(18)  is  substituted  into  eq.(19).  The  2D-strain 
compatibility  equation  is  given  by 

^^xz.xz  ~  ^z.xx  ~^x.zz  (20) 

Substitution  of  eq.(12)  into  eq.(20)  and  use  of  eqs.(17)  and  (18)  leads  to  a  fourth  order 
ordinary  differential  equation 

Kf~(z)-X^(2R[,  +  RQf;(z)+X'RiJ.(zhX%,E,  (21)- 


The  solutions  /■  (z)  can  be  expressed  by 

/,  U)  =  X  exp(m^,z, )+ 

j=\  «33 

where  Ay,  are  constant  and  the  values  of  m  ,  are  given  by 


^ m,;  ^ 


j 


=  ±A 


a-  4- b; 


V  J 


=  ±A 


a.  -  b; 


1 

V 


V  y 


where 


a  I  —  /?55  +  2/?, 3 


c,.  =  2R‘ 


The  stress  field  solutions  of  the  i'''  layer  can  be  found  using  eq.(22)  as 


(22) 


(23) 


(24) 
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It  is  noted  that  the  above  solutions  satisfy  the  boundary  conditions  of  eq.(14)  while  the 
surface  traction-free  boundary  conditions  of  eq.(15)  and  the  interface  continuity  conditions 
of  eq.(  16),  remain  to  be  satisfied,  resulting  4N  equation  for  4N  unknown  constants  Aj. . 

Figure  18  shows  some  of  the  results  obtained  by  the  2D  elasticity  model  and  FEM,  the 
distribution  of  the  inter-laminar  shear  stress  ,  and  normal  stress,  ,  near  the  free  edge 


along  the  longitudinal  (x  axis)  direction. 
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Figure  18  Distribution  of  normal  stress,  G ^  and  inter-laminar  shear  stress,  T^.  near  free  edge  of  a  standard 
bimorph  piezo-actuator  under  applied  voltage  of  208V,  predicted  by  (a)  the  2D  elasticity  model  and  (b)  FEM 
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It  follows  from  Fig.  18  that  the  predicted  interlaminar  shear  stress  near  the  free  edge  by 
the  2D  elasticity  model  is  close  to  those  by  2D  FEM  analysis. 


2.  Material  Property  Modeling  For  Fiber  Based  Piezoelectric  FGM  Bimorphs 

Our  work  on  piezoelectric  FGM  bimorphs  has  led  to  the  need  for  accurate  material 
properties  for  piezoelectric  composites.  One  possible  fabrication  method  for  piezoelectric 
FGM  actuators  is  through  the  use  of  stacked  layers  of  piezo-ceramic  fibers.  In  these 
devices  the  fiber  volume  fraction  is  graded  from  layer  to  layer  thereby  resulting  in  graded 
material  properties,  including  piezoelectric  coupling,  from  layer  to  layer.  Due  to  the  fact 
that  the  fiber  diameter  is  nearly  the  same  as  the  layer  thickness  within  the  composite  and 
that  within  an  FGM  the  fiber  distribution  is  clearly  not  uniform  through  the  thickness  of 
the  composite,  a  uniform  packing  of  fibers  in  these  devices  is  not  achieved.  Because  of 
this  deviation  from  idealized  fiber  packing  many  of  the  assumptions,  including  matrix 
connectivity,  of  most  micromechanics  models,  including  those  based  on  Eshelby’s  (1957) 
work,  are  violated.  This  leads  to  the  use  of  the  finite  element  method  (FEM)  for 
evaluation  of  individual  unit  cells.  Others  have  looked  at  this  approach,  but  they  assumed 
symmetries  which  are  not  actually  present  in  reality. 


2.1.  The  Fiber  Based  Piezoelectric  FGM  Bimorph 

The  basic  design  consists  of  stacking  two  piezoelectric  FGM’s  into  a  symmetric 
bimorph  as  shown  in  Figure  19.  The  FGM  laminates  consist  of  piezoelectric  ceramic 
fibers  contained  in  an  epoxy  matrix.  A  second  type  of  fiber,  which  is  not  piezoelectric,  is 
used  as  a  "place  holder"  in  the  layers  with  lower  piezoelectric  fiber  volume  fraction,  Vf ,  to 
prevent  adjacent  layers  from  being  pressed  into  each  other  during  the  cure  cycle. 


Poling 

Direction 


f 


•Top  Electrode  (Al) 


«m*o°.«oVrVr, 


.  Brass  Centra! 
Electrode 


^  Bottom  Electrode  (Al) 
P2T  Fibers 
Second  Fiber  Type 


Figure  19.  Schematic  end  view  of  the  fiber  based  FGM  bimorph  actuator. 
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2.2.  Modeling  Fiber  Based  Piezoelectric  Composite  Material  Properties 

For  the  case  of  current  piezoelectric  fiber  composites,  the  diameter  of  the  fibers  is 
almost  as  large  as  the  layer  thickness  and  simply  examining  Vf,  does  not  give  an  accurate 
picture  of  the  packing  of  the  fibers.  For  this  reason,  the  concept  of  ceramic  line  fractions 
in  circular  fiber  composites  was  used.  For  this  work  the  line  fractions  are  labeled  Y  and  Z, 
where  Y  is  the  ratio  of  fiber  diameter  to  fiber  centerline  spacing  within  each  layer,  and  Z  is 
the  ratio  of  fiber  diameter  to  layer  thickness. 


(A)  (B) 

Figure  20.  Rectangular  versus  square  packing  of  circular  fibers  in  a  composite.  (A)  Rectangular  packing, 
top  image  showing  Vf  =  .56,  Z  =  .99  and  Y  =  .71.  (B)  Square  packing,  top  image  showing  v,-  =  .56,  Z  = 
Y=.84. 

This  distinction  between  square  and  rectangular  fiber  packing  is  important  as  Vf  in 
the  FGM  layers  will  be  solely  a  function  of  Y,  the  in  plane  line  fraction  because  all  layers 
are  the  same  thickness.  Further,  when  considering  that  the  dielectric  constants  of  the  fiber 
and  the  matrix  differ  by  a  factor  of  100,  examination  of  the  top  images  in  Figure  20 
indicates  that  the  electric  field  will  be  very  different  in  the  two  different  packing 
arrangements. 

To  show  the  deviation  from  the  idealized  results  of  most  micromechanics  models, 
it  is  useful  to  compare  the  FEM  results  to  an  analytical  model.  The  recently  published 
work  of  Mikata  (2000)  allows  determination  of  the  Eshelby  tensor  for  a  composite  where 
the  poling  axis  is  perpendicular  to  the  long  axis  of  the  fiber.  Even  with  the  limitations 
discussed  above,  use  of  Mikata’s  Eshelby  tensor  in  the  Mori-Tanaka  (1973)  mean  field 
theory,  as  outlined  by  Dunn  and  Taya  (1993b)  will  provide  a  comparison  to  the  FEM 
models  to  be  used  for  material  property  determination. 

The  first  step  in  the  FEM  analysis  is  to  determine  which  unit  cell  to  use  and  also  which 
planes  of  symmetry  may  be  used  to  minimize  the  computational  effort  required  to  solve  the 
model.  For  this  work  the  1/2  symmetric  unit  cell  was  chosen  as  shown  in  Figure'll. 
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Figure  21.  Schematic  representation  of  1/2  symmetric  unit  cell  with  faces  labeled  for  boundary  condition 
application.  Poling  axis  is  labeled  as  P  in  the  figure.  Faces  3a  and  3b  are  the  surface  electrodes. 

The  choice  to  use  the  1/2  symmetric  unit  cell  which  models  the  entire  layer  thickness  was 
due  to  the  fact  that  there  is  not  a  plane  of  electrical  symmetry  parallel  to  the  1-2  plane 
halfway  between  faces  3a  and  3b.  This  electrical  asymmetry  is  clearly  seen  from  the  plot 
of  voltage  versus  thickness  in  a  layer  of  the  composite  as  shown  in  Figure  22. 


Figure  22.  Plot  of  Voltage  versus  thickness. 


The  symmetry  planes  in  this  model  are  faces  lb  and  2b  or  the  1-3  and  2-3  planes.  These 
are  clearly  planes  of  mechanical  and  electrical  symmetry.  The  boundary  conditions  for 
this  unit  cell  are  listed  below  in  Table  7. 
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Face 

Mechanical  B.C.’s 

Electrical  B.C.’s 

1a 

All  Ui  coupled 

D,  =  0 

2a 

All  Ug  coupled 

0 

II 

CM 

Q 

3a 

All  U3  coupled 

V  coupled 

1b 

(Symm)  u^  =  0 

p 

II 

0 

2b 

(Symm)  U2  =  0 

□2  =  0 

3b 

All  U3  coupled 

V  coupled 

Table  7.  Summary  of  boundary  conditions  on  FEM  1/2  unit  cell. 

The  model  was  meshed  using  8-node  hexahedral  elements  and  a  mapped  mesh,  as  shown 
in  Figure  23,  as  this  would  provide  more  accurate  results  than  provided  by  a  free  mesh  of 
4-node  tetrahedrons.  The  choice  of  the  hexahedral  elements  would  also  allow  easier 
control  of  mesh  density  in  areas  with  high  gradients  such  as  the  thin  layer  of  matrix 
material  between  the  fiber  and  the  surface  electrodes. 


Figure  23.  Exploded  view  of  unit  cell  showing  hexahedral  elements.  Fractions  are  Z  =  .99,  Y  =  .56,  Vf  =  .43. 
Shown  in  same  orientation  as  unit  cell  of  Figure  4. 

The  number  of  elements  used  in  this  work  ranged  from  3540  to  3660.  The  small 
range  is  due  to  the  fact  that  most  of  the  elements  are  in  the  region  directly  around  the  fiber 
and  over  the  small  range  of  Z  line  fractions  considered,  the  shape  of  the  elements  is  fairly 
constant.  The  solution  at  this  element  density  was  able  to  quickly  converge  and  there  are 
no  poorly  shaped  elements.  The  larger  number  of  elements  was  for  the  low  Y  line 
fraction.  The  small  change  added  a  few  more  elements  in  the  area  between  fibers  within 
each  layer,  although  there  was  not  too  much  concern  in  this  area  as  the  gradients  were  all 
fairly  low.  The  practical  upper  limits  of  this  model  are  Z  =  .995  and  Y  =  .97  for  Vf  =  .76, 
or  very  close  to  the  maximum  Vf  of  0.785  for  square  packed  cylindrical  fibers. 
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The  effective  material  properties  were  then  found  for  a  range  of  Y  and  Z  line 
fractions  by  using  the  FEM  method  to  run  a  series  of  hypothetical  experiments  based  on 
those  of  Bent  and  Hagood  (1997)  with  the  addition  of  several  new  cases  to  determine  the 
rest  of  the  material  properties.  These  experiments  are  summarized  in  Table  8  and  the 
material  properties  for  these  calculations  are  listed  in  Table  9.  It  is  further  noted  that  in 
cases  5  through  7,  some  of  the  symmetry  conditions  listed  in  Table  7  were  removed  and 
replaced  with  a  different  set  of  constraints  as  the  symmetric  displacements  would  not 
always  be  valid  in  the  shear  cases. 


Case 

Mechanical  Loading 

Electrical  Loading 

Properties  Obtained 

1 

Uj  =  0  on  faces  1  a,  2a,  3a,  3b  (zero  strain) 

Apply  voltage  V  on  top  electrode. 
Ground  bottom  electrode 

®31»  ®32>  ®33>  ^33 

2 

Apply  displacement  Ui  on  face  la 

Ui  =  free  on  2a,  3a,  3b 

U2  =  U3  =  0  on  1a,  2a,  3a,  3b 

V  =  0  on  Electrodes  (short  circuit) 

C11,  C12.  Ci3,  631 

3 

Apply  displacement  U2  on  face  2a 

U2  =  free  on  la,  3a,  3b 

Ui  =  U3  =  0  on  la,  2a,  3a,  3b 

V  =  0  on  Electrodes  (short  circuit) 

Ci2i  C22t  C23,  632 

4 

Apply  displacement  U3  on  face  3a 

Us  =  free  on  la,  2a;  U3  =  0on3b 

Ui  =  U2  =  0  on  la,  2a,  3a,  3b 

V  =  0  on  Electrodes  (short  circuit) 

C32»  C32,  C33,  ©33 

5 

Apply  displacement  U2  on  face  3a 

U2  =  free  on  1  a,  2a,  2b;  U2  =  0  on  3b 

Ui  =  U3  =  0  on  la,  2a,  2b,  3a 

V  =  0  on  Electrodes  (short  circuit) 

C44,  024 

6 

Apply  displacement  U3  on  face  la 

U3  =  free  on  2a,  3a,  3b;  U3  =  0  on  1  b 

Ui  =  U2  =  0  on  la,  1b,  2a,  3a,  3b 

V  =  0  on  Electrodes  (short  circuit) 

^55»  ®15 

7 

Apply  displacement  on  face  2a 
u,  =  free  on  1a,  1b,  3a,  3b;  U3  -  0  on  2b 

U2  =  U3  =  0  on  la,  1b,  2a,  2b,  3a,  3b 

V  =  0  on  Electrodes  (short  circuit) 

^66 

Table  8.  Load  cases  for  FEM  1/2  symmetric  unit  cell. 


C11 

(Gpa) 

C12 

(Gpa) 

Ci3 

(Gpa) 

C33 

(Gpa) 

C44 

(Gpa) 

®31 

(C/m^) 

®33 

(C/m^) 

®1S 

(C/m^) 

£11/80 

£33/^0 

PZT-5A 

121 

75.4 

75.2 

111 

21.1 

-5.4 

15.8 

12.3 

916 

830 

Improved 

Epoxy 

10.8 

5.7 

5.6 

10.9 

2.7- 

-0.002 

0.007 

0.005 

7.3 

7.3 

Table  9.  Material  properties. 


Some  of  the  results  of  these  analyses  are  plotted  below  in  Figures  24  and  25  along  with  the 
results  from  the  micromechanics  model  using  Mikata’s  Eshelby  Tensor  as  mentioned 
above.  The  dashed  vertical  line  is  Vf  =  0.785,  or  the  maximum  square  packing  of  round 
fibers. 
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PZT  Fiber  Volume  Fraction 


0.2  0.4  0.6  Q.i 

PZT  Fiber  Volume  Fraction 


Figure  24,  Effective  material  properties  (Cn,  C12,  C33,  C55,  631,  633)  versus  Vf. 
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Figure  25.  Effective  £33/60  versus  Vf.  (A)  Compares  FEM  results  to  epoxy  (white)  matrix  Eshelby.  (B) 
Compares  FEM  results  to  PZT  (black)  matrix  Eshelby  (Note:  reversal  of  volume  fractions  here.) 

From  examination  of  Figures  24  and  25  it  is  seen  that  the  elastic  stiffness  constants 
were  independent  of  different  Z  line  fractions  and  were  predicted  reasonably  well  by  the 
micromechanics  model.  However,  for  the  piezoelectric  and  dielectric  constants  a  strong 
dependence  on  the  Z  line  fraction  was  observed,  but  it  was  also  noted  that  the  results 
trended  towards  the  micromechanics  model  as  the  Z  line  fraction  was  reduced.  Based  on 
this  strong  dependence  of  the  piezoelectric  and  dielectric  constants  on  the  Z  line  fraction, 
it  is  important  to  minimize  the  thickness  of  the  matrix  layer  between  the  fibers  and  the 
electrode. 

Comparison  of  Figures  25(A)  and  25(B)  gives  insight  into  the  "black  matrix/white 
matrix"  issue  over  which  species  is  actually  the  matrix  at  high  Vf.  The  white  matrix  case 
gives  the  general  lower  limit  for  all  but  the  highest  Z  line  fraction  and  the  black  matrix 
case  gives  an  upper  limit  for  Z  =  1 . 

2.3.  Analysis  Of  The  Piezoelectric  Fiber  Based  FGM  Bimorph 

Once  the  effective  composite  material  properties  had  been  found  it  was  then 
possible  to  predict  the  response  of  a  fiber  based  FGM  bimorph.  From  the  analyzed  cases, 
the  middle  Z  line  fraction  case  of  Z  =  .99  was  chosen  at  a  Y  line  fractions  of  Y  =  .93,  .71, 
.50  to  represent  layers  which  are  all  PZT  fibers,  2/3  PZT  fibers  and  1/2  PZT  fibers 
respectively.  The  actuators  as  modeled  then  had  three  100  |im  thick  piezoelectric 
composite  layers  on  each  side  of  a  25  pm  thick  brass  central  electrode  witfi  10  pm  thick 
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aluminum  electrodes  on  the  surface  for  a  total  thickness  of  645  |im  and  a  total  length  of  50 
mm.  This  represents  the  type  of  actuators  that  are  currently  being  fabricated.  These  were 
modeled  using  ANSYS  in  the  2D  plane  stress  mode  which  gives  displacements  that  agree 
well  with  a  full  3D  model,  but  with  much  less  computational  effort.  The  results  for  this 
FGM  bimorph  are  compared  to  a  non-FGM  bimorph  which  was  dimensionally  identical, 
but  had  the  effective  properties  for  the  Z  =  .99,  Y  =  .93  layer  throughout.  Both  of  these 
models  were  subjected  to  an  applied  electric  field  of  500  volts/mm.  These  results  for  axial 
stress,  Qx,  are  shown  below  in  Figure  26. 

From  these  results  it  is  seen  that  the  FGM  bimorph  is  superior  with  a  tip 
displacement  187  |im,  and  a  rnaximum  axial  stress  in  the  piezoelectric  layers  of  1.27  MPa 
compared  to  175  pm  and  2.08  MPa  respectively  in  the  non-FGM  bimorph.  The  reduced 
stresses  are  likely  due  to  the  gradual  build  up  in  stress  through  the  three  layers  of  the  FGM 
instead  of  the  much  quicker  ramp  up  evidenced  in  the  single  layer  of  the  non-FGM 
samples.  The  FGM  bimorph  also  compares  very  favorably  with  the  results  for  a 
monolithic  FGM  monomorph  reported  elsewhere,  as  the  tip  displacement  of  our  FGM 
bimorph  at  500  volts/mm  equates  to  37.4  pm/cm  bending  displacement  which  is  superior 
to  the  values  of  28  pm/cm  for  the  monolithic  FGM  monomorph. 

Initial  comparison  of  the  FGM  to  non-FGM  bimorphs  may  only  show  a  7% 
increase  in  tip  deflection,  but  there  are  other  benefits  to  the  FGM  which  go  along 
favorably  with  this  increase.  Use  of  the  FGM  bimorph  as  opposed  to  the  non-FGM 
bimorph  has  led  to  a  reduction  in  axial  stress  of  39%.  In  addition  to  this,  there  is  the  fact 
that  the  FGM  only  uses  about  75%  of  the  expensive  ceramic  piezoelectric  fibers  required 
for  a  non-FGM  bimorph  of  identical  thickness,  resulting  in  significant  cost  savings  and 
reduced  weight. 


(A) 


(B) 


Figure  26.  Results  for  axial  stress,  (A)  FGM  bimorph.  (B)  non-FGM  bimorph. 


Further  examination  of  the  plots  in  Figure  26  show  that  the  highest  stresses  are 
actually  in  the  aluminum  surface  electrodes.  With  this  in  mind,  it  will  be  possible  to 
achieve  higher  displacements  by  replacing  the  aluminum  foil  with  a  different  electrode 
material  that  has  a  lower  stiffness  (more  flexible),  such  as  a  conductive  polymer  film  or 
conductive  epoxy  paint. 

3.  Design  of  Laser  Based  Displacement  Measurement  System 

To  allow  accurate  measurement  of  tip  displacement  induced  in  the  FGM  bimorph  a 
laser  based  displacement  measurement  system  was  designed  and  fabricated  as  shown 
below  in  Figure  27.  The  use  of  a  laser  based  transducer  allows  for  non-contact 
measurement  of  tip  displacement  to  allow  for  true  no  load  displacement  rheasurements. 
The  incorporation  of  and  x-y  translation  stage  into  the  measurement  system  allows 
measurements  to  be  easily  taken  at  various  points  on  the  specimen.  The  data  acquisition 
and  control  of  the  system  is  done  with  a  Windows  based  PC.  The  accuracy  of  the  system 
is  ±3|xm  in  one  location  and  ±5jj,m  when  scanning  over  a  surface  using  the  x-y  stage. 
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Figure  27.  UW  Designed  and  fabricated  laser  based  non-contact  displacement  measurement  system. 


4.  Concluding  Remarks  And  Future  Work 

The  use  of  FGM  in  piezoelectric  laminates  reduces  the  stress  field  while  increasing 
the  out-of  plane  displacement  as  compared  with  the  standard  bimorph.  The  optimization  of 
the  material  property  gradient  profile  within  the  FGM  laminae  provides  higher  out-of¬ 
plane  displacement  while  reducing  the  stress  field.  The  CLT  and  2D  elasticity  solution  of 
out-of-plane  displacement  and  stress  field  are  close  to  those  obtained  by  FEM  analysis. 

The  effective  material  properties  for  piezoelectric  fiber  based  composites  were 
predicted  by  using  an  FEM  model  of  the  fiber  unit  cell  which  allowed  rectangular  fiber 
packing  to  be  accounted  for.  It  was  found  that  the  elastic  constants  are  mostly  only  a 
function  of  the  overall  volume  fraction  Vf,  while  the  piezoelectric  and  dielectric  constants 
were  more  strongly  affected  by  changing  the  Z  line  fraction  alone  than  by  changes  in  Vf. 

Using  these  material  properties  allowed  modeling  of  piezoelectric  fiber  based 
FGM  bimorph  actuators  which  are  under  development  and  indicates  that  they  have 
advantages  over  non-FGM  bimorphs  and  also  have  higher  deflection  than  for  other 
existing  monolithic  (non-fiber  based)  FGM  monomorphs. 

Fabrication  of  such  actuators  is  considered  to  be  a  key  target  of  our  future  work,  as 
well  as  the  validation  of  the  analytical  modeling  as  compared  to  the  experimental  data.  In 
conclusion,  the  concept  of  FGM  is  found  to  be  very  useful  in  reducing  the  stress  field  at 
the  interface  where  large  mismatch  in  strain  exists,  thus  inducing  large  stress  field,  which 
can  be  reduced  by  using  the  FGM  interface  design. 
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